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ABSTRACT 


Water-type  solar  heating  systems  have  received  much 
attention  in  the  literature  and  several  good  computer  simulation 
models  are  available.   Air-type  space  heating  units  have 
advantages  which  may  make  them  more  suited  to  a  particular 
need  than  the  water-type  system. 

This  computer  simulation  offers  a  model  of  an  air-type 
solar  heating  system  with  which  one  can  determine  the  monthly 
system  output.   All  significant  system  parameters  are  easily 
varied  as  is  the  storage  bed  size  and  material.   The  simulation 
accepts  all  air-type  collector  designs.   The  residence  is 
simply  modeled  and  estimating  its  heat  load  is  reduced  to 
calculating  the  product  of  floor  area  times  a  constant. 
Insolation  data  input  is  based  on  monthly  averaging-techniques. 

The  quantities  varied  in  the  present  application  of  the 
model  were  collector  size  and  storage  volume. 
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NOMENCLATURE 


Numbers  in  (.  )  refer  to  the  equation  (,s)  in  which 
the  symbol  appears . 


A  (12,14,15)        cross-sectional  area  of  pebble  bed  storage 

unit     m2 

A  overall  collector  area     m2 

c 

A   (13)  total  surface  area  of  pebble  bed  storage 

p  unit     m2 

A  „  (9)  area  of  collector  normal  to  air  flow     m2 

xs 

b  subscript  refering  to  the  storage  bed 

material 

c,   (13,15)  specific  heat  of  pebbles  in  storage  bed; 

assumed  constant  at  0.837  KJ/Kg  -  °C 

cf  (11-13,15)       specific  heat  of  air;  assumed  constant 

at  1.012  KJ/Kg  -  °C 

D  (14)  diameter  of  pebbles,  all  assumed  equal   m 

D„  (7,9)  characteristic  length  of  collector 

(hydraulic  diameter)  for  use  in  determining 
Reynold's  number;  equal  to  twice  the 
spacing  between  the  plates     m 

d  (23)  dust  factor 

F'  (11)  collector  efficiency  factor  [1] 

FD  (10,11)  collector  heat  removal  factor  [1] 

f  subscript  refering  to  the  fluid  (air) 

G  (14)  flow  rate  per  unit  cross-sectional  area 

of  storage  bed     Kg/s  -  m2 

G   (11)  flow  rate  per  unit  of  collector  area 


c 


Kg/s  -  m 


2 


H  (16,2  3)  monthly  average  daily  total  radiation  on  a 

horizontal  surface     KJ/m2  -  day 
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h,  (7)  convective  heat  transfer  coefficient  for 

collector  absorber  plate     W/m2  -  °C 

h0  (7)  convective  heat  transfer  coefficient  for 


l2 


air 


kins 

(5) 

L 

(12 

,15) 

m 

(12 

,13) 

m 

(9, 

12,13,15) 

N 

CI, 

25)  : 

(20) 

Nu 

I    (7 

,8) 

n 

(12 

,13, 

15) 

nc 

,    (27) 

collector  back  plate     W/m  -  °C 


h   (6)  radiation  heat  transfer  coefficient  between 

the  two  air  duct  surfaces  of  collector 
(sides  neglected)     W/m2  -  °C 

h   (12,14,15)       volumetric  heat  transfer  coefficient  for 

storage  bed  material     W/m3  -  °C 

K,  (16,23)  ratio  of  the  monthly  average  daily  diffuse 

radiation  to  the  daily  extra-terrestrial 
radiation 

K.  (16,23)  ratio  of  the  monthly  average  daily  total 

radiation  to  the  daily  extraterrestrial 
radiation 


(7)  thermal  conductivity  of  air,  taken  to  be 

constant  at  0.029  W/m  -  °C 

thermal  conductivity  of  insolation  used 
under  collector,  taken  to  be  constant  at 
0.043  W/m  -  °C 

length  of  storage  bed     m 

subscript  referring  to  the  mth  segment 
of  the  storage  bed 

mass  flow  rate  of  air     Kg/s 

number  of  glass  covers:  day  of  year  measured 
from  January  1 

Nusselt  Number 

number  of  segments  into  which  the  storage 
bed  is  considered  to  be  divided  for 
application  of  numerical  technique 

refractive  index  of  glass 
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Qu  (10)  hourly  rate  of  useful  thermal  energy 

collection  per  unit  area  of  collector 
surface     W/ra2 

R.   (16,17,23)       ratio  of  radiation  on  a  tilted  surface 

to  that  on  a  horizontal  surface 

R      (8,9)         Reynold's  Number  for  collector  air  flow 
D,H  based  on  hydraulic  diameter 

r,  (16,21,23)       ratio  of  hourly  to  daily  diffuse  radiation 

r.  ratio  of  hourly  to  daily  total  radiation 

S  (19,23)  hourly  flux  of  solar  radiation  absorbed 

by  collection  surface     W/m2 

SI  (16)  hourly  flux  of  solar  radiation  striking 

a  tilted  surface     W/m2 

S  (23)  shading  factor 

T  ,   (1,10)         ambient  temperance     °C 

T   .  (10)  inlet  air  temperature  to  the  collector    °C 

c ,  1 

TV   (12,13)         temperature  of  pebbles  contained  in  storage 
m  segment    m 

Tf   (12,13)         temperature  of  air  entering  storage  segment 
m  m     °C 

Tf     (12,13)       temperature  of  air  leaving  storage  segment 
m+11  m+1     °C 

Tr  mean  fluid  temperature  in  collector 

r  ,m 

T   (1)  mean  absorber  plate  temperature  during 

"  operation  of  collector     °K 

T   (13)  temperature  of  the  surroundings  containing 


s 


the  storage  bed    °C 


T  (6)  mean  temperature  for  radiation  between 

the  absorber  plate  and  back  plate  of 
collector,  assumed  constant  at  340 °K 
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T,  (6)  mean  temperature  of  absorber  plate    °K 

T2  (6)  mean  temperature  of  back  plate     °K 

tins  ^  thickness  of  insulation  behind  back  plate  m 

U,   (5)  overall  collector  back  loss  coefficient 

b  W/m2  -  °C 

UT  (10,11)  overall  collector  loss  coefficient 

L  W/m2  -  °C 

UT   (13)  overall  storage  loss  coefficient    W/m2  -  °C 

La 

U. (45°)  (1,4)       overall  collector  top  loss  coefficient 

for  a  unit  tilted  45   from  the  horizontal 
W/m2  -  °C 

U.  (3)  (4)  overall  collector  top  loss  coefficient 

for  a  unit  tilted   degrees  from  the 
horizontal     W/m2  -  °C 

V  (3)  wind  velocity,  assumed  constant  at  5  m/s 
a  solar  altitude 

a  (6  )  (24,28)       directional  absorptivity  of  collector  plate 
3  (4,16,23)         collector  tilt  from  horizontal 

Y  (18)  surface  azimuthal  angle  (due  soth  being 

zero  and  east  positive) 

6  (18-20,22)  declination 

e  (13,15)  void  ratio  of  storage  bed 

£  (1)  emissivity  of  glass  cover  plates 

£  (1)  emissivity  of  absorber  plate 

XT 

e,  (6)  emissivity  of  absorber  plate 

£2  (6)  emissivity  of  back  plate 

0.  (17,18,24-27)     angle  of  incidence  of  direct  radiation 

measured  from  the  normal  to  the  collector 
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6  (26,27)  angle  of  refraction 

y  (9)  absolute  viscosity     Kg/m-s 

p  (13,15)  density  of  pebble  bed  material     Kg/m3 

p    (16,23)  diffuse  ground  reflectivity 

p(6  )  (25,26)       reflectivity  of  a  single  air-glass  inter- 
face for  direct  radiation  as  a  function 
of  4>t 

p (60  )  (24)         reflectivity  of  a  single  air-glass 

interface  for  diffuse  radiation 

time     s 

f.e  .  1  effective  transmittance  -  absorptivity 

e   '  >    (23)      product  for  direct  (beam)  radiation 

and  diffuse  radiation,  respectively 

t (0. )  (24)  directional  transmittance  of  glass  as  a 

function  of  0. 

4>  (18,19,22)        latitude  (north  positive) 

oj  (18,19,21)        hour  angle  (solar  noon  zero  and  morning 

positive) 

oj   (21,2  2)  hour  angle  corresponding  to  sunrise 

(and  sunset  since  the  day  is  assumed 
symmetric  about  solar  noon) 


a  =  5.6693(10)  8  W/m2  -  °KI+ 
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I.   INTRODUCTION 

This  study  was  undertaken  for  the  purpose  of  developing 
a  computer  simulation  of  an  air-type  residential  solar 
space-heating  model  from  which  to  determine  the  optimum 
collector  and  storage  size  for  a  given  heat  load.   Furthermore, 
the  need  for  comparison  between  this  and  liquid-type  systems 
dictated  that  the  model  developed  be  adaptable  to  changes 
in  various  significant  parameters  and  components. 

To  date  the  majority  of  published  computer  simulation 
models,  as  well  as  actual  installations,  has  been  of  the 
"liquid-type",  that  is,  a  solar  heating  system  wherein 
water  or  a  water-antifreeze  solution  is  the  working  fluid  [7]. 
The  most  noteable  exception  to  this  trend  has  been  the 
simulations  and  installations  of  Dr.  G.  Lof. 

Liquid-type  systems  are  favored  [7]  for  commercial 
installations  where  a  large  amount  of  energy  must  be 
distributed.   In  such  a  case  the  storage  tank,  piping  and 
insulation  required  occupies  far  less  valuable  space  than 
would  the  storage  and  ductwork  necessary  to  transport  the 
same  amount  of  heat  using  air  at  reasonable  rates  of  flow. 
Also,  the  basic  components  (heat  exchangers,  pumps  and 
controls)   associated  with  the  liquid-type  system  become 
proportionally  less  expensive  as  the  size  of  the  system 
they  serve  grows . 
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On  the  scale  of  a  single- family  dwelling  however,  it 

appears  a  trade-off  becomes  practicable  between  liquid  and 
air  systems  depending  upon  the  geographic  location,  nature 
of  an  existing  conventional  system  (for  retrofit  solar 
installations)  and  extent  of  solar  energy  utilization  desired, 
among  other  considerations.   Some  salient  trade-off 
considerations  of  the  two  types  of  systems  currently  most 
widely  employed  are  summarized  in  Table  1.   Though  not 
exhaustive,  this  indicates  generally  some  of  the  more 
significant  differences  between  the  systems. 

These  points  serve  only  as  a  first-cut  at  making  an 
overall  system  evaluation  because  various  design  strategies 
may  be  employed  to  lessen  the  disadvantages  of  the  liquid 
system  or  extend  the  utilization  of  the  air  system. 

For  example,  if  in  a  water-type  system  the  space 
heat  is  supplied  by  blowing  air  over  a  heat  exchanger 
directly  into  the  room, instead  of  relying  on  natural 
convection  from  baseboards  or  radiators,  lower  water 
temperatures  may  be  acceptable.   Solar  system  efficiency 
increases  significantly  when  the  thermal  energy  is  absorbed 
at  the  lowest  useful  temperature.   This  heat  transfer 
can  be  accomplished  at  the  storage  container  (e.g.,  in  a 
rock  bed  or  finned  annulus  surrounding  the  tank)   with 
standard  ducting  distributing  the  heated  air  throughout  the 
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TABLE  li   Comparison  of  relative  merits  of  linuid  and  air-type 
■olar  heating  systems.   "Liquid"  refers  water  or 
water-antifreeze  solution  as  the  thermal  energy 
transport  medium  coupled  to  a  water  storage  tank. 
"Air"  refers  to  air  as  the  transport  medium  with  a 
rock  bed  for  storage. 

The  base  for  comparison  is  space  heating  with 
additional  abilities  treated  as  advantages  over  this. 


LIQUID  AIR 


ADVANTAGES 

1.  smaller  storage  volume  per 
unit  of  heat  stored 

2.  transport  network  occupies 
smaller  space 

3.  directly  adaptable  for  use 
with  absorption  A/C 
equipment 

4.  directly  adaptable  to 
domestic  hot  water  supply 

5.  adaptable  to  either  forced- 
hot-water  or  forced-hot-air 

6.  if  chillocl- water  air-condi- 
tioning is  supplied,  one 
distribution  system  serves 
for  both  heating  and  cooling 


DISADVANTAGES 

more  conplex  controls,  heat 
exchangers  required  in  most 
applications,  maintainance 
more  costly 

Less  efficient  energy  collection 
at  temperatures  required  for 
forced-hot-water  baseboard 
heating  (140-190OF) 

may  require  some  treatment  of 
transport  liquid,  such  as: 
antifreeze,  antiscale, 
antieorrosive  or  biocide 


possibi lity 
in  dwelling 
leakage,  or 
freeze-up 

water  may  be 
antifreeze  ( 
uncertain  fu 
availability 

where  centra 
is  supplied, 
heat  require 
system  (pipe 


of  expensive  damage 

due  to  collector 
to  collector  duo  to 

scarce  commodity; 
oil  by-product)  has 
ture  cost  and 

1  air  conditioning 

forced- hot- water 
s  parallel  distribution 
s)  for  space  heating 


ADVANTAGES 

1.  simple  operation,  longer  lifetime, 
less  maintainance,  less  capita), 
and  installation  cost 

2.  since  the  transport  fluid  is 
identical  to  the  ultimate  delivery 
medium,  no  heat  exchanger,  per  :e, 
is  necessary 

3.  system  can  be  adapted,  though  not 
readily,  to  cooling  of  dwelling* 

4.  where  central  air  conditioning  is 
supplied,  one  distribution  system 
serves  for  both  heating  and  cooling 
(ductwork  sized  for  ;_he  ni&aS  flwW 
rates  associated  with  A/C  lends 
itself  to  collection  of  heat  at  the 
lowest  temperature) 

5.  air  is  free,  rocks  are  plentiful 
and  leaks  are  neither  dangerous 
nor  damaging 


DISADVANTAGES 

not  readily  adaptable  (possibly  not 
feasible)  to  providing  source  of 
heat  for  conventional  air 
conditioning  equipment 

where  (finned  tube)  heat  exchangers 
are  used,  they  must  be  larger  than 
that  required  in  a  liquid-system 
to  transfer  the  same  quantity  af: 
energy 

ducting  and  storage  reouire  more 
space  for  the  same  energy  transport 
or  storage,  respectively 


'e.g.,  cool  storage  at  night  for  use 
during  the  day. 
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dwelling.   Alternatively  the  heated  water  can  be  piped  to 
those  locations  where  space  heat  is  required  and  air  blown 
over  a  coil- type  heat  exchanger  into  the  space.   The  latter 
scheme  might  be  advantageous  where  space  is  at  a  premium 
or  where  the  distribution  network  is  extensive. 

On  the  other  hand,  by  placing  a  water  coil  within  the 
rock  bed  of  the  air-type  system  or  girdling  the  bed  with 
such  a  coil,  preheated  or  fully  heated  domestic  hot  water 
can  be  provided.   Alternatively,  by  cooling  the  rock  bed 
by  night  with  ambient  air,  the  dwelling  can  be  cooled 
during  the  day  using  the  storage  as  a  heat  sink.   The  materials 
of  the  rock  bed  being  rather  low  in  cost,  it  might  be 
worthwhile  to  provide  a  separate,  uninsulated  storage  bed 
a  meter  or  so  underground  to  be  cooled  at  night  in  the 
above  manner  while  also  losing   heat  to  the  ground;  this 
unit  would  be  exclusively  for  summer  air  conditioning 
and  sized  accordingly.   The  main  storage  unit,  well  insulated 
and  connected  to  the  collector,  would  be  used  to  heat 
domestic  hot  water  (as  described  above)  and  would  be 
available,  as  well,  to  provide  space  heating  as  needed 
on  cool  summer  nights  or  days.   This  latter  bed  would  be 
sized  for  winter  space  heating. 

As  can  be  seen,  hybridization  and  innovation  can  be 
used  to  great  advantage,  especially  in  adapting  to  an 
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already  existing   conventional  system,  but  makes  a  clear-cut 
determination  of  the  "best"  system  for  a  given  situation 
very  difficult.   Regional  weather  characteristics  and  labor 
rates,  the  cost  in  dollars  of  one  system  versus  another 
(including  various  levels  of  component  refinement  in  each 
case)  and  the  availability  of  materials  (e.g.  a  reliable 
source  of  water  or  antifreeze)  only  complicate  further  a 
fundamentally  complicated  problem  in  thermodynamics.   To 
this  end,  "canned"  computer  simulations  appear  to  offer  the 
most  direct  and  cost-effective  manner  of  gleaning  from 
the  myriad  of  interrelationships  the  optimum  design  for  a 
given  dwelling  and  location. 
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II.   DESCRIPTION  OF  THE  SOLAR  HEATING 
SYSTEM  SIMULATED 

A  schematic  diagram  of  the  solar  space  heating  system 
employed  in  this  study  appears  in  Figure  1.   The  residence 
is  heated  by  forced-hot-air.   The  flat  plate  collector 
transforms  the  absorbed  portion  of  the  incident  solar 
radiation  into  thermal  energy  which  is  then  extracted  by  the 
working  fluid  passing  over  the  heated  absorbing  plate 
(see  Figure  2) .   This  heat  transfer  is  accomplished  with 
air,  the  ultimate  delivery  medium.   When  satisfying  immediate 
heating  requirements  of  the  dwelling,  air  circulates 
directly  between  the  collector  and  dwelling  -  cycle  A  -  if 
useful  heat  is  available.   When  no  solar  energy  is  directly 
available,  the  room  air  is  circulated  through  the  storage 
medium  -  cycle  B  -  if  useful  heat  is  available  there.  Should 
there  be  a  requirement  for  space  heating  beyond  the  capacity 
of  both  these  processes,  auxiliary  heat  -  cycle  C  -  is 
provided  to  make  up  the  deficit.   Finally,  when  solar 
energy  is  available  for  which  there  is  no  immediate  need, 
the  working  fluid  is  circulated  through  the  storage  medium 
transferring  thermal  energy  thereto.   As  structured,  the 
absence  of  independent  heat  exchangers  (the  pebble  bed  is 
at  once  the  storage  medium  and  a  single  heat  exchanger) 


on 
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insulcition 


C  auxiliary 
source 


storage 
loss 


residence 


CYCLE 
A 
B 


D 


heating  dwelling  with  collector 

heating  dwelling  with  storage  (no 
useful  heat  available  from  collector 
directly) 

heating  dwelling  with  auxiliary  source 
(not  enough  useful  heat  available  from 
collector  or  storage) 

Heating  storage  with  the  collector 
(dwelling  requires  no  heat) 


FIGURE  1:   Schematic  diagram  of  the  solar  space 
heating  system  used  in  this  study. 
Operational  modes  are  listed  in  order 
of  precedence. 
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disallows  the  ability  to  charge  and  discharge  storage 
simultaneously  as  is  possible  in  liquid-type  systems.   Since 
the  transport  fluid  can  also  be  used  as  the  delivery  medium, 
however,  no  advantage  is  lost. 

Auxiliary  heat,  necessary  for  those  periods  when  no 
direct  or  stored  solar  energy  is  available  from  the  system, 
is  provided  by  a  conventional  heater.   Mounted  in  the  outlet 
duct  as  shown  in  Figure  1,  this  heat  exchanger  (electrical 
resistance  or  gas/oil  fired)  has  an  auxiliary  inlet  in  order 
that  circulation  through  storage  does  not  occur.   An 
exception  to  this  logic  could  be  the  case  where  off-peak 
electrical  rates  make  charging  storage  with  auxiliary, 
off-peak  energy  economical  thus  supplementing  storage  input 
from  the  collector. 

The  storage  bed  is  considered  to  be  located  within 
the  heated  space,  therefore  convective  heat  losses  from 
the  unit  are  regarded  as  an  uncontrolled  gain  to  the  room. 

As  shown,  the  location  of  the  circulation  fan  would  cause 
a  reduced-pressure  in  the  collector.   Collectors  for  such 
systems  are  difficult  to  seal  and  this  configuration  would 
allow  infiltration  of  cold  air.   However,  this  has  the 
effect  of  lowering  the  overall  collector  temperature 
resulting  in  more  efficient  operation  at  the  expense  of 
lowering  the  collector  outlet  temperature.   If  one  intends 
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to  collect  the  thermal  energy  at  the  lowest  useful 
temperature  anyway,  this  method  would  appear  to  offer  an 
advantage.   This  simulation  does  not  involve  such 
considerations,  nor  does  it  account  for  pressure-drops 
through  storage  or  elsewhere.   In  this  regard,  it  is 
assumed  peripheral  calculations  have  addressed  compatability 
between  mass  flow  rates  used  and  various  flow  areas  in  the 
system  (collector  and  ducting  cross-sections,  pebble-bed 
effective  flow  cross-section,  etc.). 
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III.   SIMULATION  OF  THF  SYSTEM  COMPONENTS 

A.    Collector 

Several  flat-plate  collectors  are  available  for  air-type 
solar  heating  systems.   This  simulation  is  written  to 
accept  such  units  when  the  collector  characteristics  have 
been  reduced  to  expressions  for  tL  and  F' .   The  derivation 
of  these  terms  and  other  factors  for  a  given  collector 
geometry  is  given  in  [1].   The  particular  collector 
geometry  chosen  for  this  simulation  is  shown  schematically 
in  Figure  2. 

To  reduce  computing  costs,  the  empirical  relationship 
developed  by  Klein  for  the  top  loss  coefficient,  U.  (45°) 
was  used: 


Ut  (45°)  = 


N  1 

+  — 


<344/Tp)[(Tp-Tamb)/(N+f)]°-31    1^ 


-1 


g<yTamb)(Tp2+Tamb2> 

[e  +0.0425N(l-e  ) I"1  +  [ (2N+f-l) /e  ]-N 
P  p  y 


(1) 


where 


f  =  (1.0  -  0.04  hw  +  5.0  (10 )"4  hw2  )  (1  +  0.058  N)     (2) 
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i°f 


V,T. 


glass  covers  (2) 


flow 


absorber  plate 
back  plate 


J£L 


T 


T 


U, 


el/hl,Tl 


e2,h2'T2 


F'    = 


Ut   +   Ub 


U, 


1   + 


h,    + 


1  T   l/h2  +   l/hj 


FIGURE  2:   Schematic  of  Collector  Used  in  This 
Simulation  (Elevation) 
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and 

h  =  5.7  +  3.8  V  (3) 

(See  NOMENCLATURE  section  for  definitions) 

As  this  expression  corresponds  to  a  collector  tilt 
angle,  6,  of  45 °,  an  additional  equation  to  interpolate 
for  any  tilt  angle  (and,  incidentally,  for  plate  emissivities 
other  than  0.95)  was  employed  [1] : 


Ut(B)  =  Ut(45°)  [1  -  (3-45)  (0.00259  -  0.00144)  e  )]   (4) 

where  f  is  dimensionless,  V  is  in  m/s ,  all  temperatures 
are  in  °K  and  hT,  U.  (45°)  and  U.  (3)  are  in  units  of 
W/m2-°C.       V   t  •  t 


In  equation  (1)  the  mean  plate  temperature,  T  ,  is  not 
known.   As  a  result,  an  iterative  solution  is  used  in  the 
determination  of  U,  (45°),  (see  Appendix).    Wind  velocity,  V, 
is  taken  as  a  constant  year  round  but  is  easily  varied  if 
desired.   The  expression  used  for  the  back  loss  coefficient 
is: 

Thus  the  overall  collector  loss  coefficient,  UL,  is 
determined  to  be: 


UL  =  Ut  +  Ub  (6) 
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The  term  hr,  representing  the  radiative  heat  transfer 
coefficient  between  the  back  and  absorbing  plates  (the 
sides  of  the  duct  being  neglected  in  the  heat  transfer) ,  is 
determined  from  the  righthand-most  expression  below: 

o(T1+T9) (T,2+T,2)       4aT3 

=  (7) 

l/e1+l/e2-l      l/e1+l/e2-l 


where  T,  the  mean  temperature  for  radiative  heat  transfer,  is 
taken  as  340 °K-  This  assumption  is  not  critical  for  the 
varying  conditions  encountered  in  normal  operation  [1]. 


The  convective  heat  transfer  terms  -  h-,  and  h2  -,    taken 
to  be  equal,  are  based  on  a  fully  developed  forced  turbulent  air  flow 
between  the  parallel,  flat  plates.  For  the  mass  flow  rate  and  all 
the  collector  sizes  used  in  this  simulation,  the  Reynold's 
Number  (based  on  the  hydraulic  diameter)  was  greater  than 
2100  and  the  length-to-diameter  ratio  is  larger  than  10, 
thereby  justifying  the  assumption.   If  this  is  not  the  case, 
one  must  consider  the  effect  of  the  developing  boundary 
layer  (see  Section  4.13  of  [1]).   The  resultant  expression 
is: 

h]_  =  h2  =  Nu  (*air/DH)  ,   W/m2-°C  (8) 
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where 


Nu  =  0.0158 


m  D 


H 


Axs  » 


0.8 


(9) 


With  these  relations  and  that  for  IL,  F*  is  determined 


F'  = 


U, 


1  + 


hl  +  l/h0  +  1/h, 


(10) 


l2    -'  ~r 

Depending  upon  the  current  temperature  of  the  room  or 
storage  "cool  side"  (either  of  which  might  be  the  collector 
inlet  temperature  depending  on  the  operating  mode  in  effect) 
the  energy  absorbed  by  the  collector  may  or  may  not  be 
available  at  a  useful  temperature.   The  expression  for  useful 
energy  gain  from  the  collector  per  unit  area  in  either  case 
is  [1]  i 


Qu  -  FR  IS  "  VTc,i  "  Tamb>l 


(11) 


where   F^  is    the   collector  heat  removal   factor  equal   to    [1] : 


G  C 

Fp   =   &    II   -   exp(-    UTF'/G   C_)J 

R  UT  P 


(12) 


In  the  simulation,  if  the  derived  quantity  Qu  is  negative 
or  zero,  the  collector  is  not  operated.   Assuming  S,  the 
absorbed  solar  energy  (see  equ.  (24))  is  positive  the  inter- 
pretation of  a  negative  value  of  Qu  is  that  if  run,  the  collector 
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would  be  operating  at  an  average  temperature  too  high 
relative  to  the  energy  being  absorbed  to  result  in  an 
increase  in  fluid  temperature.   Expressed  differently,  more 
thermal  energy  would  be  lost  to  the  surroundings  than  is 
being  absorbed  by  the  collector. 

Due  to  the  long-term  nature  of  the  energy  balance 
simulated,  the  thermal  capacitance  of  the  collector  was 
neglected.   This  is  justified  because  the  thermal  energy 
"lost"  in  bringing  the  collector  to  operating  temperature 
in  the  morning,  is  subsequently  returned  to  the  system  at 
the  end  of  the  solar  collection  period.   If  the  time 
response  of  the  system  were  of  interest,  for  control  purposes 
for  example,  this  simplification  would  not  be  valid. 

Neglected  also  was  a  consideration  of  the  trade-off 
involved  between  collecting  marginal  amounts  of  useful 
energy  and  the  electrical  energy  expended  by  the  circulation 
fan  to  do  so. 

Collector  tilt  was  taken  to  be  about  20°  greater  than 
the  latitude  which  emphasizes  that  this  simulation  intends 
to  size  the  collector  for  midwinter  heating.   Between  the 
practical  range  of  tilt   angles    (latitude  +  23.45°), 
particular  angles  favor  particular  collection  schemes. 
For  example,  $  =  (latitude  +  23.45°)  results  in  the  largest 
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ratio  of  absorbed-to- incident  solar  energy  occuring  on 
December  21  -  this  favors  space  heating.   If  3  =  (latitude  - 
23.45°) ,  the  largest  ratio  occurs  on  June  21  and  favors 
energy  collection  on  that  day. 

This  simulation  assumes  an  absorber  plate  emissivity 
for  long  wave  radiation  of  0.95.   With  selective  surfaces 
available  at  higher  cost,  this  value  can  be  depressed  while 
maintaining  reasonable  solar  absorptance  values.   When 
mass-production  techniques  are  employed  to  reduce  the  cost, 
the  significant  increase  in  absorbed  energy  should  justify 
the  additional  expense. 

The  effect  of  dust  on  the  coverplates  is  accounted 
for  as  is  shading  due  to  the  collector  frame  work.   The 
number  of  covers  is  taken  to  be  two.   Ground  reflectance  [1] 
is  taken  as  0.2  year-round  which  is  conservative  because 
a  value  of  0.7  is  suggested  when  there  is  snow  cover. 

See  Section  III. E .  for  the  equations  used  to  determine 
S,   the  absorbed  portion  of  the  radiant  energy  SI  incident 
upon  the  outermost  coverplate  of  the  titled  collector. 
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B.    Storage 

The  pebble  bed  storage  model  is  taken  from  Mumma  and 
Marvin  [2J.   Two  difference  equations  comprise  the  mathematical 
model  with  one  derived  from  an  air-side  energy  balance 
and  the  other  from  the  pebble-side.   The  equations  used 
are  [2]  : 


Tf     =  (Tf  -Tb  )exp[-(hv  AL)/m  Cf  n) ]  +  Tb 
m+1      mm  m 


(13) 


and 

T   (x+Ax)  =  [iCf(Tf»(T)%l(T)1"nLA\(T)-TB)1 
bm  (pAL/n)  (l-e)Cb 


+  T,  (t) 
m 


(14) 


where  subscript  f  refers  to  the  air,  subscript  b  refers  to  bed 
material  (rock) ,  t  is  time  and  subscript  m  refers  to  one  of 
the  n  segments  into  which  the  storage  bed  is  divided  for  using 
the  model.   Other  terms  are  defined  in  NOMENCLATURE. 


The  term  h   is  the  volumetric  heat  transfer  coefficient 


defined  by  [1] : 


hv  =  650 


G/A 


0.7 


W 

~3 

m  -°C 


(15) 
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The  procedure  to  use  these  equations  is  [2J: 

(a)  Assume  an  initial  bed  temperature  distribution 
at  time  equal  zero. 

(b)  For  a  specified  bed  inlet  fluid  temperature, 
evaluate  the  fluid  temperature  entering  and 
leaving  all  n  sections  of  the  bed  using 
equation  (13 ) . 

(c)  Evaluate  the  new  bed  temperatures  at  each 
section  after  a  time  increment  At  using 
equation  (14 ) . 

(d)  Step  forward  in  time  by  At. 

(e)  Return  to  step  (b)  and  repeat. 

The  "specified  bed  inlet  fluid  temperature"  is  either 
the  room  temperature  or  collector  outlet  temperature 
depending  upon  the  operating  mode.   The  time  step,  At  in 
seconds,  can  be  any  value  between  0  and  3600  in  this 
simulation,  corresponding  to  the  various  operational 
profiles  which  are  possible.   In  order  to  insure  that  the 
solution  be  stable,  the  time  step  in  use  must  satisfy 
the  following,  [2]: 


(pAL/n) (l-e)c. 

At  <   - S. (IS) 

m  cf  [l-exp(-hvAL/m  df  n)] 


on 
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The  computer  simulation  used  automatically  determines 
the  largest  number  of  storage  sections,  n,  which  will 
insure  a  stable  solution  for  the  largest  possible  time 
step  which  may  occur  -  3600  sec. 

In  a  well-designed  storage  bed  the  pebble  size  is  small 
enough  that  the  temperature  gradient  in  the  individual  pebbles 
is  insignificant  [1].   Physically  this  means  the  only 
significant  resistance  to  heat  transfer  to  or  from  the 
storage  medium  is  that  due  to  the  convective  heat  transfer 
phenomenon  occuring  at  the  surface  of  the  pebbles. 

The  convective  heat  transfer  coefficient  associated 
with  the  surface  of  the  pebbles  represents  a  thermal 
resistance  over  which  we  have  little  control.   On  the  other 
hand,  we  can  manipulate  to  some  extent  the  thermal  response 
of  the  interior  of  the  pebbles  by  advantageously  choosing 
the  size  and  shape  of  the  pebble  in  relation  to  its  material 
composition.   It  is  generally  accepted  [3]  that  if  the  Biot 
number  (representing  the  ratio  of  internal  resistance  to 
external  heat  transfer  resistance)  is  less  than  0.10/ then 
computational  error  resulting  from  assuming  a  uniform 
temperature  within  the  pebble  will  be  less  than  5%.   This 
assumption  has  been  made  here  in  order  to  represent  the 
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volumetric  heat  transfer  coefficient,  h  ,  by  equation  (15) . 
If  the  pebble  size  is  to  be  varied  to  examine  the  effect, 
this  should  be  taken  into  account. 

The  void  ratio,  e,  is  chosen  to  be  0.3.   This  corres- 
ponds to  a  body-centered-cubic  packing  arrangement  of 
exactly  similar  spheres.   Practically,  this  is  the  smallest 
ratio  which  might  occur. 

Early  trial  simulations  included  unsteady  conduction  in 
order  to  destratify  the  temperatures  in  the  storage  segments 
while  storage  was  idle.   The  standard  explicit  numerical 
technique  used  13]  introduced  unacceptable  errors  due  to 
the  relatively  large  longitudinal  dimension  of  the  storage 
segments.   This  length  is  dictated  by  other  stability  require- 
ments.  An  implicit  method  may  have  improved  the  results  but, 
no  such  capability  is  presently  included  in  this  simulation. 
Though  less  elegant,  it  is  believed  no  misrepresentation  of 
the  storage  response  occurs  due  to  neglecting  this  effect. 

During  the  heating  season,  taken  here  to  be  September 
through  May  inclusive,  the  storage  bed  is  assumed  to  be 
located  within  the  residence  such  that  convective  losses 
from  storage  are  not  "lost"  at  all  but  contribute  to  the 
thermal  energy  requirements  of  the  dwelling. 
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For  the  months  of  June,  July  and  August  it  is  assumed 
these  losses  take  place  to  the  ambient  temperature. 
In  actual  practice  the  storage  bed  would  probably  be 
located  in  the  basement  of  the  dwelling  where  some  of  the 
convective  loss  would  be  recovered  but  a  significant 
portion  would  be  lost  to  the  ground.   In  order  to  avoid 
unnecessary  complication,  the  rather  artificial  device  of 
locating  the  storage  bed  based  on  the  season  was  used  in 
order  not  to  unreasonably  misrepresent  one  operating 
mode  more  than  the  other. 

All  the  fluid  and  storage  material  (rock)  properties 
were  considered  constant  at  mean  system  operating 
temperature values  (see  Table  2) . 
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C.    Residence 

The  residence  is  modeled  as  a  degree-day  heat  sink. 
It  was  chosen  to  correspond  in  that  characteristic,  and  size, 
to  that  used  by  Butz  [4]  in  his  water-type  system  (17,O00Btu/ 
°F-day  =  374  w-hr/°C-hr) .   Assumed  to  be  a  box-like  structure, 
it  has  dimensions  of  45'  x  40'  x  8'  with  standard 
construction  materials,  3  inches  of  insulation,  2  x  4's 
on  16-inch  centers  and  15%  fenestration.   The  effects  of 
heat  loss  through  the  floor,  infiltration  of  outside  air, 
solar  radiation  through  unshaded  windows  and  internal  heat 
generation  from  occupants  and  electrical  appliances  were 
considered  by  Butz.   Considering  only  the  walls,  roof  and 
windows  as  media  through  which  conductive  and  convective 
heat  transfer  may  take  place,  one  arrives  at  a  figure  of 
18,200.0  Btu/ °F-day;  a  7%  error  which  overpredicts  the  energy 
need.   It  appears  a  much  less  strenuous  calculation  may  be 
adequate  to  characterize  a  residence. 

The  model  produces  data  based  on  the  following 
assumptions  regarding  a  desirable  home  heating  scenario: 

(i)   During  the  heating  season  the  comfort  range  is 
between  67.0-70.0 *F  (19.4  -  21.1  °C)  between  the  hours 
of   0630  and  2230  and  62.0  -  65.0°F  (16.7  -  18.3  °C)  for 
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the  remainder  of  the  24  hour  period.   For  the  purpose  of 
determining  storage  loss  to  the  surroundings  or  collector 
inlet  temperature  during  direct  residence  heating  from  the 
collector,  the  mean  of  the  desired  temperature  range  in 
effect  is  used.   It  is  tacitly  assumed  that  the  residence 
is  within  the  desired  temperature  range  at  all  times  by 
some  combination  of  direct,  stored  or  auxiliary  energy. 

(ii)  During  June,  July  and  August  the  desired  room 
temperature  is  between  70.0  -  73.0°F  (21.1  -  22.8  °C) 
between  hours  of  0630  and  2230  and  65.0  -  73.0°F  (18.3  - 
22.8°C)  for  the  remainder  of  the  24  hour  period.   Although 
during  these  months  storage  temperatures  are  allowed  to 
be  180 °F  (a  temperature  too  high  to  be  acceptable  for  space 
heating),  space  heating  would  in  fact  be  possible  with  a 
simple  mixing  damper  in  the  ductwork.   The  effect  would  be 
to  divert  only  part  of  the  fixed  air  flow  through  storage. 
The  present  model  does  not  account  for  this  possibility  and 
simply  provides  auxiliary  heat  when  necessary  during  these 
summer  months. 

(iiL)  If  the  ambient  temperature  is  greater  than  the 
desired  room  temperature  but  less  than  the  maximum  acceptable 
room  temperature  (73.0°F)  no  input  of  any  kind  is  provided. 
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D.    Auxiliary  Heat 

Auxiliary  heat  is  supplied  in  the  amount  necessary 
to  make  up  the  deficit  between  that  which  the  solar  heating 
system  is  able  to  provide  (including  uncontrolled  storage 
losses  to  the  room)  and  that  required  by  the  residence. 
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E*    Insolation  and  Weather 

The  method  of  Liu  and  Jordan  [5]  was  used  in  this 
simulation  to  provide  insolation  input.   The  development 
and  application  of  this  method  follows  directly  that  of 
Felske  [6].   In  this  approach,  actual  hourly  values  of 
insolation  for  an  entire  year  are  unnecessary.   Instead, 
this  method  requires  only  that  the  monthly  average  values 
of  H   Krp  and  K,  be  determined  for  a  given  location.   These 
terms  are  defined  in  NOMENCLATURE  and  [5]  explains  how 
they  are  derived.   The  location  used  in  this  study  was 
arbitrarily  chosen  to  be  New  York  City  (40.77°^). 

With  the  above  parameters  determined  for  each  month, 
one  is  able  to  calculate  for  an  average  day  of  that  month, 
the  hourly  flux  of  solar  radiation  incident  upon  the 
collector  from  the  equation  [6J: 


SI  =  rdH{[Rb(l-Kd/Kt)]  +  [1/2  (l+cos3) Kd/Kt] 

V ^- J  ' v ' 

beam  diffuse  sky 


+  £1/2  (l-cos3)p   ]}  (17) 

v H *_• 

diffuse  ground  reflection 


where  it  is  assumed  that  r.  =  r.  and: 
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cos    6, 

Rb      =     (18) 

sin  a 


cosG.    =   cos6    costo    (cosy   sin<{)   sinB   +  cost))   cosft) 
+   siny    sing   cos6    sinoo 
+   sin6  (sin<{>   cosg   -   cosy   cost})    sin$)  (19) 

sina  =   sin6    sin<(>   +   cos6   cosif)   cosu  (20) 

<5    =    23.45   sin    [360    (284   +  N)/365J  (21) 


IT  COSU)     -     COSOJ 

rd  = £_  (22) 

24    sinw„    -    w   cosoo 
s  s  s 


coso3_   =  -   tan<{>   tanfi  (23) 


where  N  is  here  the  number  of  consecutive  days  to  mid-month, 
measured  from  January  1.   When  the  quantity  SI  is  negative 
or  zero  the  collector  was  not  operated. 


Of  course   the  insolation  striking  the  tilted,  outermost 
surface  of  the  collector  is  not  the  energy  which  ultimately 
reaches  the  absorbing  plate.   The  hourly  rate  of  radiation 
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striking  and  being  absorbed  by  the  collector  plate  (accounting 
for  dust,  shading,  reflections  and  absorption)  was  computed 
from  the  relation  [6]: 


s  =  (Ta)eff,b  ra*l%   U-V*t>J 


(24) 


rdH 


+  (xa)eff^d  -^-[ (l+cose)Kd/Kt  +  (l-cos8)p  J  (1-s)  (1-d) 


where  from  [6] : 


(xa) 


x(6t)  a(0t) 


eff 


l-[l-a(0t)]p(6O°) 


(25) 


and 


T(ftt)  = 


l-p(6t) 


l+(2N-l)p(6t) 


(26) 


P(6t)  =  2 


sin2  (62"et)    tan2  (62"et) 


sin2  (62+6t)    tan2  (e2+et) 


(27) 


-1 

■  sin 


sine 

t 

_Ng 

(28) 


in  which  N  equals  the  number  of  glass  cover  plates  and  other 
terms  are  defined  in  NOMENCLATURE. 
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Reference  [6]  indicates  that  the  directional  absorptivity 
of  the  collector  surface,  a(0.),  can  be  represented  by  a 
curve  for  a  typical  surface.   The  equation  provided  by  Felske 
to  represent  this  typical  curve  is: 

20.03 

a(G.)  =  1.0 r-T  (29) 

(100.03  -  et) 

where  6.  assumes  its  actual  hourly  value  for  beam  radiation 
but  is  assigned  the  value  of  60°  when  used  to  represent  the 
diffuse  radiation  component  [6]. 

The  hourly  ambient  temperature  distribution  is  more 
random  than  insolation  and  it  can  probably  be  argued  that 
one  reasonable  distribution  of  the  average  daily 
temperature  is  about  as  good  as  another.   In  this  simulation, 
the  scheme  illustrated  in  Figure  3  was  used.   The  average 
daytime  temperature  and  nighttime  temperature  were  cast 
on  a  sine  wave  peaking  at  1300  hours  (Note:   throughout  this 
simulation,  hourly  values  are  referenced  to  solar  time) . 
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FIGURE  3:  Distribution  of  Ambient 
Temperature  Through  the 
Day 
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TABLE  2 
Simulation  Design  Parameters 


Insulation  thickness,  collector 

k  . 
air 

kins 
m 

D 

N 

ng 

Plate  spacing,  collector 

s 

V 

a(60°) 

3 
e 


0. 

837  KJ/Kg-°C 

1. 

012  KJ/Kg-°C 

0. 

02 

0. 

076m 

0. 

029  W/m-°C 

0. 

043  W/m-°C 

1. 

2  Kg/s 

0. 

0508m 

2 

(number  of  covers) 

1. 

,526 

0. 

,  01m 

0, 

.03 

5, 

.  0  m/s 

0, 

.8346 

60° 

0, 

.3 

0 

.88 

0 

.95 

0 

.0 
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y  (absolute  viscosity) 

P 

V 


1.912  CIO)"5  Kg-m/s 

40.77°N 

2400  Kg/m3 

0.2 
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TABLE  3 
Approximate  Time  Constants  for  Major  System  Components 

Values  for  collector  and  storage  are  based  on  the  basic 
component  parameters,  assumed  constant  at  approximately  the 
mean  operating  temperature  associated  with  each  component. 

The  value  for  the  residence  is  calculated  based  on 
the  time  to  change  the  interior  temperature  of  the  author's 
home  a  given  number  of  degrees  over  a  period  of  time  with 
constant  ambient  temperature  and  no  heat  input. 

Collector:         1/4  hour 
Storage  Bed:       1  hour 
Residence:         10  hours 
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IV.   OPERATIONAL  LOGIC  AND  STRATEGY 

The  flow  chart  of  Figure  4  represents  the  operation  of  the 
solar  system.   A  breakdown  of  the  computer  logic  used  to 
implement  this  basic  logic  appears  in  the  Appendix. 

This  simulation  operates  by  accounting  for  the  following 
heat  transfers  from  and  into  the  residence  unit: 

a.  The  heat  loss  from  dwelling  driven  by  the 
difference  between  the  desired  room  temp- 
erature and  ambient  temperature. 

b.  Heat  supplied  directly  to  the  dwelling  from 
the  collector  on  an  as-needed  basis  when 
available. 

c.  Heat  supplied  to  the  dwelling  from  the  storage 
bed  in  the  form  of  either  uncontrolled  con- 
vective  losses  or  forced-hot-air  plus  the 
uncontrolled  losses. 

d.  Heat  supplied  from  the  collector  to  the 
storage  bed. 

e.  Auxiliary  heat  as  necessary  to  meet  what 
remains  of  the  heat  load  should  sources 
(b)  and  (c)  prove  insufficient  to  meet  the 
entire  demand. 
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NOTE: 


Energy  cannot  be 
supplied  to  room 
by  both  collector 
and  storage 
during  the  same 
hour. 


(COLLECTOR  AREAS 
STORAOf,  VPI.UMry 


FIGURE  4  Flow  Diagram  of  Overall  Simulation  Logic 
(Refer  also  to  Figure  A. 2) 
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f.    Total  useful  thermal  energy  obtaining  from 
the  collector  for  all  purposes. 

The  only  temperatures  monitored  or  used  as  controlling 
input  are  those  of  the  storage  bed  segments.   The  temperature 
of  the  residence  is  assumed  at  all  times  to  be  within  the 
comfort  range  in  effect  for  a  given  season. 

To  implement  the  insolation  representation  of  Liu  and 
Jordan  [  5 ] ,  it  has  been  assumed  that  the  monthly  average 
of  the  hourly  insolation  and  ambient  temperature  applied  for 
one  24-hour  period  predicts  the  system  behavior  for  the 
entire  month  when  the  results  of  this  "average"  day  are 
multiplied  by  the  number  of  days  in  the  month.   Butz  [  4 ] 
used  a  similar  extrapolation  with  hourly  Weather  Bureau  data. 

A  discontinuity  occurs  in  the  system  operation  in  going 
from  one  month  into  another  due  to  the  change  in  these  monthly 
average  values.   Furthermore,  the  final  temperature  profile 
in  the  storage  bed  carries  into  the  succeeding  month.   For 
these  reasons  the  first  days  of  operation  in  any  given  month 
do  not  represent  a  true  average  daily  response.   If  the 
simulation  is  operated  for  several  days  utilizing  the 
insolation  and  ambient  temperature  profiles  of  that  month, 
the  influence  of  the  starting  conditions  is  eliminated.   A 
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measure  of  the  degree  to  which  the  starting  conditions  have 
been  obscured  is  the  correspondence  between  storage  bed 
temperature  profiles  separated  by  24  hours  of  operation. 

For  the  ranges  of  parameters  and  values  of  constants 
used  here,  4  to  5  days  operation  served  to  damp  out  starting 
effects  in  most  cases.   Therefore,  the  following  strategy 
was  used: 

i.    For  each  month,  9  full  days  of  hourly 

operation  were  simulated, 
ii.   The  last  run  (corresponding  to  the  9th 
day)  is  used  as  the  "average"  day.   The 
results  obtained  for  this  day  are  multiplied 
by  the  number  of  days  remaining  in  the  month, 
iii.  The  results  for  the  entire  month  are  obtained 
by  adding  the  cumulative  results  of  the  first 
9  days  to  those  obtained  in  (ii) . 

During  the  heating  season,  the  maximum  allowable  average 
temperature  of  the  storage  bed  is  chosen  to  be  145°F.   In 
the  limit,  then,  this  is  the  highest  possible  forced-air 
temperature  which  can  be  derived  from  storage.   This  temp- 
erature is  at  least  20°F  below  that  which  outlets  from  a 
conventional  forced-hot-air  heating  unit  where  common 
practice  indicates  a  rise  in  temperature  of  100 °F  or  more 
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through  the  furnace  18] .   This  lower  operating  temperature 
allows  more  efficient  operation  of  the  collector  but  requires 
a  higher-than-conventional  mass  flow  rate  of  warm  air  into 
the  room.   The  codes  18]  restrict  the  rate  of  room  air  change- 
out  and  must  be  consulted.   If  upon  supplying  the  hourly 
room  needs  the  collector  is  able  to  bring  storage  to  this 
temperature,  the  collector  ceases  to  operate.   Subsequent 
collector-to-storage  operation  is  only  to  the  extent  necessary 
to  maintain  this  temperature  in  the  face  of  demands  made  upon 
storage.   This  mode  of  operation  causes  the  low  usefulness 
ratios  (useful  heat  collected/incident  thermal  energy)  in 
early  fall  and  late  spring. 

During  the  summer  months  when  this  simulation  does  not 
supply  solar  space  heating,  the  maximum  storage  temperature 
is  allowed  to  be  180 °F  for  the  purpose  of  domestic  water 
heating.   Higher  temperatures  are  possible  if  pressure  relief 
and  mixing  values  are  used  in  the  water  system  but  then  the 
collector  would  operate  less  efficiently.   180°F  was  chosen 

to  avoid  unnecessary  complication.   No  attempt  is  made  to 
model  a  realistic  daily  hot  water  usage  scenario  in  the  summer. 
Instead,  the  heating  load  presented  to  the  bed  is  that  of  an 
increased  convective  heat  loss  coefficient  corresponding 
approximately  to  an  80  gal/day  draw-down  of  100 °F  water  spread 
evenly  through  the  24-hour  day, added  to  the  previously  used 
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coefficient.   Accurate  assessment  of  the  water  heating  mode 
does  not  contribute  to  the  objective  of  this  study;  this 
portion  of  the  simulation  is  included  only  for  continuity. 

If  the  collector  can  heat  the  dwelling  directly  but  is 
unable  to  meet  the  full  requirement,  auxiliary  heat  is  used 
to  make  up  the  deficit.   In  this  situation  the  only  contri- 
bution from  storage  is  the  uncontrolled  convective  loss. 
This  strategy  is  used  in  order  to  obtain  all  the  useful 
solar  energy  available  and  does  not  consider  the  electrical 
energy  required  by  the  circulation  fan  to  obtain  this 
marginal  amount  of  thermal  energy. 

Several  arbitrary  choices  have  been  made  as  to  what 
constitutes  a  desirable  system  operation. 

i.    If  when  heating  storage  with  the  collector  the 
increase  in  average  bed  temperature  corresponds 
to  a  thermal  energy  input  of  less  than  3.3  W-hr 
in  one  minute  (an  hourly  rate  of  200  W)  ,  the 
system  is  not  operated. 

ii.   When  heating  the  room  from  storage  a  1%  overshoot 
of  the  room's  needs  is  tolerated.   The  overshoot 
shows  up  as  excess  heat  which  is  subtracted  from 
the  otherwise  useful  energy  gain  of  the  system 
and  results  from  the  time  step  used  (see  iv) . 
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iii.  If  at  any  time  while  heating  the  room  from 

storage  the  thermal  energy  added  in  one  minute 
is  less  than  5%  of  l/60th  the  total  hourly  room 
need,  then  the  system  is  relaxed  and  the 
auxiliary  source  provides  the  remainder  of 
the  heat  load. 

iv.   All  heat  transfers  to  and  from  storage  are 

accomplished  in  1-minute  intervals.   If  some 
fraction  of  a  minute  over  the  hour  is  used  or 
an  unacceptably  large  overshoot  in  energy  input 
occurs,  the  1-minute  interval  in  which  this 
occurred  is  nullified  and  the  process  continues 
in  5- second  intervals.   Any  further  overshoot 
resulting  during  such  a  reduced  interval  is 
accepted  and  appears  as  an  accountable  excess. 

When  heating  storage  from  the  collector,  the  inlet 
temperature  to  the  collector  is  held  fixed  for  the  duration 
of  the  hour  at  its  value  at  the  beginning  of  the  hour. 
In  fact  this  temperature  increases  as  the  bed  is  heated 
and,  if  strictly  accounted  for,  could  result  in  a  decrease 
in  useable  thermal  energy  available  to  storage.   Errors  due 
to  this  assumption  are  considered  small  considering  the 
small  change  in  temperature  of  the  storage  bed  which  occurs 
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in  one  hour.   So,  in  the  interest  of  decreasing  computation 
time/  the  effect  was  neglected. 
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V.   RESULTS 

Figure  5  demonstrates  the  ability  of  various  collector 
sizes  to  meet  the  heating  demands  of  the  residence  used 
in  this  simulation.   As  collector  size  increases  a  greater 
proportion  of  the  room  needs  can  be  accommodated.   If 
analysis  of  cost  tradeoffs  and  regional  weather  character- 
istics indicate  that  it  is  advantageous  to  use  solar  power 

for  approximately  half  the  residential  heating  needs  in 

2 

mid-winter,  then  this  figure  indicates  that  the  60  m 

collector  should  be  employed. 

The  curves  of  the  ratio  of  solar  energy  obtained  to 
room  needs  for  various  collector  sizes  were  derived  using 
a  fixed  storage  volume.   Figure  6  indicates,  however,  that 
mid-winter  system  capability  for  various  collectors  is 
relatively  insensitive  to  the  size  of  storage.   The  reason 
for  this  result  is  that  for  the  collector  sizes  examined, 
no  collector  provides  on  a  daily  basis  more  energy  than 
the  residence  consumes  in  that  time.   There  is  then  no 
net  accumulation  of  energy  in  storage  and  the  total  of  the 
daily  solar  energy  obtained  reflects  the  room  needs  only. 

The  uppermost,  flat  curve  of  Figure  5  represents 
the  practical  limit  of  system  output.  Achieving  this 
output  requires  a  collector  of  such  size  that  storage  is 
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brought  to  a  fully  charged  state  each  day  while  the 
immediate  daytime  dwelling  heat  load  is  met  as  well.   A 
"fully  charged"  storage  bed  is  one  at  the  maximum  allowable 
average  temperature.   It  is  assumed  in  presenting  this 
limiting  case  that  sufficient  storage  volume  has  been 
assigned  to  the  system  to  meet,  in  full,  the  overnight 
heat  load  if  the  bed  is  initially  fully  charged. 

In  Figure  5  the  QU/QROOM  curve  (see  SIMULATION 
NOMENCLATURE  for  definition  of  terms)  for  the  largest 
collector  used  is  seen  to  behave  differently  than  the 
other  two  curves  during  April  and  May.   The  "cross  over" 
observed  results  when  the  storage  bed  approaches  the 
maximum  average  temperature  allowed.   When  this  occurs, 
very  little  further  solar  energy  input  to  the  system  takes 
place. 

A  QU/QROOM  ratio  exceeding  the  maximum  practical 
value  of  1  is  seen  to  occur  in  late  spring  for  all  three 
collector  sizes  examined.   This  is  because  in  those  months 
enough  solar  energy  can  be  obtained  from  the  collector  to 
meet  the  room  needs  and  allow  a  net  accumulation  of  energy 
in  storage.   Note  that  as  collector  size  increases,  the 
QU/QROOM  curve  approaches  the  limiting  curve. 
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The  curves  of  QU/QI  of  Figures  5  and  6  represent  the 
ratio  of  solar  energy  collected  to  that  incident  upon  the 
tilted  collector  surface.   This  ratio  is  a  measure  of  collec- 
tion efficiency.   For  fixed  system  mass  flow  rate  of  air  and 
storage  size,  increasing  collector  area  decreases  collection 
efficiency  because  of  the  higher  mean  temperature  of  collec- 
tion and  storage  (even  if  panels  are  added  in  parallel) . 

Efficiency  is  generally  lower  in  early  fall  and  late 
spring  because  the  higher  insolation  and  storage  temperature 
approaching  saturation  cause  a  higher  mean  temperature  of 
collection.   This,  in  turn,  drives  a  higher  loss  rate  to  the 
surroundings  despite  relatively  warm  ambient  temperatures. 
In  late  fall  and  early  spring  the  relation  between  ambient 
temperature,  available  insolation  and  the  lower  average  temp- 
erature in  storage  combine  advantageously  to  produce  the 
highest  efficiency.   In  mid-winter  the  low  ambient  temper- 
ature drives  large  collector  losses  despite  the  rather  low 
mean  temperature  of  collection  (resulting  from  the  lowest 
insolation  and  storage  temperature  levels  of  the  year) . 

In  Figure  6 ,  the  higher  efficiency  of  larger  storage  beds 
in  spring  reflects  the  increasing  level  of  collector-to-storage 
heating  activity  and  lower  collection  temperatures  which  result 
from  distributing  about  the  same  amount  of  collected  energy 
throughout  a  larger  mass  of  rock.  The  approach  saturation  which 
will  occur  as  storage  size  decreases,  rapidly  increases  the 
efficiency  disparity. 
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Figures  7  and  8  illustrate  the  monthly  totals  of 
hourly  solar  energy  supplied  by  the  system  in  excess  of 
that  required.   During  early  spring  and  late  fall  the 
excess  is  primarily  caused  by  the  uncontrolled  convective 
storage  loss.   During  the  winter  when  all  the  storage  loss 
is  of  use  to  the  dwelling,  any  excess  occuring  is  caused 
by  an  "overshoot"  of  storage  energy  input  to  the  room 
during  active  heating  with  storage.   This  overshoot  results 
from  the  time  increment  used  in  the  simulation  (see  pg.  50 
Sec.  IV)  and  is  not  a  characteristic  of  solar  systems. 

The  excess  in  spring  and  fall  is  a  real  effect  if 
storage  is  contained  within  the  dwelling.   These  losses 
are  annoying  and  are  counterproductive  to  any  air- 
conditioning  employed  during  the  day. 

All  three  curves  of  Figures  .7  and  .8  extend  to  about 

5 
1.4(10)   W-HR  in  September,  off  the  scale  presented.   This 

occurs  because  storage  was  taken  to  be  fully  charged  at 

the  commencement  of  this  month.   The  high  convective  heat 

losses  which  are  a  consequence  cause  the  high  excess 

energy  value. 

Such  excesses  do  not  occur  in  May  because  storage 

temperature  levels  do  not  naturally  reach  saturation 

with  the  available  solar  energy  and  room  requirements 
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are  higher  in  May  than  September.   Starting  the  month  of 
September  with  a  fully  charged  storage  bed  simulates  having 
used  storage  during  the  summer  for  domestic  water  heating 
at  the  maximum  allowable  average  bed  temperature  of  180 °F 

The  two  larger  collectors  are  seen  to  provide  similar 
amounts  of  energy  in  excess  of  that  needed  for  the  month  of 
May.   This  results  because,  again,  the  storage  bed  is 
approaching  the  saturation  level.   Once  near  or  at  this 
level,  a  larger  collector  will  not  result  in  more  storage 
losses. 

In  Figure  8,  an  apparent  anomoly  occurs  in  late  spring 
when  the  excess  energy  for  the  largest  storage  unit  falls 
below  that  of  the  other  two  units.   An  examination  of  the 
hour-by-hour  energy  transfers  to  the  room  from  storage 
indicates  that  the  "overshoot"  flaw  is  the  cause  (see  pg.  50, 
Sec.  iv) .   A  lesser  percentage  of  excess  energy  input  to 
the  dwelling  occurs  each  hour  with  the  largest  bed  compared 
to  the  two  smaller  units.   This  behavior  is  evidence  of  the 
lower  overall  storage  temperature  as  bed  size  increases,  all 
else  being  equal. 

The  arbitrary  choice  of  storage  size  used  in  this 
simulation  happens  to  span  the  range  over  which  an  inversion 
occurs  in  excess  energy  output  from  storage.   It  is  observed 
in  Figure  8  that  from  October  through  December  the  excess 
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energy  first  increases  with  increasing  storage  volume  then 
decreases.   Again,  inspection  of  hour-by-hour  energy 
transfers  indicates  the  cause.   The  smallest  bed  can  provide, 
in  the  first  several  hours  after  sundown,  more  hourly  energy 
output  than  larger  units  which  have  a  lower  average  temper- 
ature.  This  is  so  despite  the  fact  that  all  beds  contain 
about  the  same  total  stored  energy. 

The  higher  hourly  output  from  the  smallest  storage  bed 
results  in  more  hourly  overshoot.   But  since  this  bed  is 
depleted  sooner,  it  is  less  often  able  to  provide  heat.   The 
cumulative  overshoot  is  consequently  less  than  for  larger 
units. 

The  largest  unit  operates  to  heat  the  room  the  most 
number  of  times,  but  each  time  the  overshoot  is  small  because 
of  the  relatively  lower  temperature  of  the  energy  provided. 
Cumulatively,  less  total  overshoot  occurs  than  with  the 
mid-size  bed  but  more  than  that  occurring  with  the  smallest. 

Figure  9  shows  the  degree  in  which  room  energy  needs  are 
met  with  the  desired  collector  and  various  storage  volumes. 
This  figure  applies  to  the  New  York  area  for  September  and 
October  starting  with  a  storage  bed  at  room  temperature. 
For  the  collector  size  and  residence  of  this  example,  the 
only  months  during  which  more  solar  energy  is  provided  than 
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that  required  to  meet  the  daily  room  needs  is  September  and 
May.   Either  of  these  months  may  be  used  to  "sense"  the 
presence  of  the  storage  bed,  no  longer  masked  by  the  over- 
powering effect  of  room  needs.   Relative  stroage  size  plays 
very  little  role  during  mid-winter  because  then  storage  only 
serves  to  distribute  solar  energy  stored  during  the  day 
through  part  of  the  immediately  following  night  with  no 
carry-over  of  energy  into  the  following  day.   It  is  clear 
that  some  storage  volume  is  necessary  but  that  size  has  yet 
to  be  determined. 

The  dashed,  vertical  line  of  Figure  9  indicates  the 
collector-to-storage  ratio  for  which  no  more  solar  energy  is 
collected  daily  by  the  system  than  is  required  for  an 
average  day  in  September.   Allowing  nothing  for  "cloudy 
days"  carry-over  storage  capability,  this  system  should  have 
a  storage  volume  in  September  equal  to  the  collector  area 
times  the  factor  1/35  (the  dimensions  being  in  meters) . 

A  storage  sizing  criterion  which  optimizes  mid-winter 
use  of  the  solar  heating  system  would  seem  to  be  a  more 
rational  criterion  than  the  above.   It  has  been  shown, 
however,  that  mid-winter  system  operation  is  relatively 
insensitive  to  storage  size.   Within  the  range  of  storage 
sizes  examined,  the  ratio  of  storage  output  to  room  needs 
increases  less  than  1%  (see  Table  4)  when  storage  volume  increases 
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twofold.   So,  such  a  criterion  is  elusive  at  best.   Intuiti- 
tively,  there  are  limits  to  the  range  of  useful  storage 
temperatures  and  this  suggests  an  upper  and  lower  bound  to 
storage  volume  exists. 

If  an  upper  temperature  limit  is  placed  on  the  bed 
(as  dictated  by  the  codes  as  regards  delivery  temperature 
or  room  air  change-out  rate,  etc. :   in  this  simulation 
145°  F  was  used)  there  is  a  lower   storage  volume  which 
will  attain,  on  average,  this  saturation  temperature  using 
all   the  energy  available  to  storage  on  a  typical  mid-winter 
day.   However,  because  the  collector  operates  inefficiently 
due  to  the  high  temperature  it  receives  from  storage,  the 
total  energy  stored  is  less  than  could  have  been  obtained 
for  any  larger  storage  bed  given  the  same  conditions. 

The  lower  limit  on  storage  bed  temperature  is  that 
which  provides  a  comfortable  influx  of  heating  air  and 
depends  on  relative  humidity.   Generally,  normal  body 
temperature  represents  a  good  estimate  of  this  lower  limit 
on  storage  temperature.   The  storage  volume  for  which  this 
temperature  obtains  when  all  the  energy  available  to  storage 
on  a  mid-winter  day  is  accepted,  is  the  upper  limit  on 
storage  volume.   This  provides  the  most  efficient  system 
operation  based  on  collected  energy  but  the  energy  stored 
at  this  temperature  is  only  marginally  useful. 
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Lacking  a  definitive  raid-winter  storage-sizing  method, 
it  is  suggested  that  the  " September  storage  sizing  criterion" 
above  be  adopted  as  a  method  of  determining  storage  volume 
given  a  collector  area. 

In  this  simulation,  using  a  5m  x  12m  collector  the 

3 

storage  volume  thus  indicated  is  1.7m  .   If  it  is  assumed 

that  the  energy  available  to  storage  in  January  for  a 

2 
6.fcm  x  2.5m  storage  volume  (see  Table  A.l)   remains 

3 
constant  as  the  storage  volume  decreases  to  1.7m  ,  the 

resulting  average  storage  bed  temperature  attains  115  F 
(if  it  starts  at  room  temperature) .   Of  course,  due  to 
collector  inefficiency  the  collected  energy  would  be  less 
as  would  be  the  temperature  attained. 

Thus,  an  even  smaller  storage  bed  could  be  used  with- 
out saturation  occurring  but  less  total  energy  would  be 
stored. 
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TABLE  4 

Energy  input  to  room  from  storage,  between  sunrise 
and  sunset,  as  a  percentage  of  room  needs  during 
that  period.   Values  are  for  January  with 
combinations  of  storage  and  collector  sizes  as 
indicated. 
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VI.   CONCLUSIONS 

From  Figures  5  and  6  it  can  be  concluded  that  collector 
area  requirements  to  meet  a  known  energy  need  can  be  deter- 
mined without  regard  to  storage  size. 

It  is  unnecessary  -co  consider  increasing  storage 
size  to  allow  for  the  "carry  over"  of  solar  energy  into 
subsequent  "cloudy  days"  because  for  economically  sized 
collector  areas  designed  for  mid-winter  needs  there  is, 
on  the  average,  no  energy  available  for  subsequent  days. 
Only  if  more  than  100%  of  the  daily  room  requirements 
can  be  supplied  from  the  system  in  mid-winter  is  there, 
on  the  average,  a  "carry  over"  capability  which  can  be 
enhanced  by  increasing  storage  size.   Collector  size  is 
the  single  most  important  consideration  in  installing 
this  solar  heating  system. 

Once  the  collector  area  necessary  to  meet  the  energy 
need  has  been  determined  a  compatible  storage  volume  must 
be  derived.   "Compatible  storage  volume"  is  taken  to  be 
that  volume  which  results  in  a  QU/QROOM  ratio  of  1  in 
September  (see  Figure  9) .   It  must  be  realized  that  such 
compatability  of  collector  and  storage  applies  only  to 
a  system  providing  approximately  50%  of  the  mid-winter 
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space  heating  requirement  and  the  location  for  which  this 
analysis  was  carried  out,  that  is,  New  York. 

Within  the  above  conditions,  it  is  suggested  that  the 
ratio  35:1  be  used  to  relate  collector  area  to  storage 
volume  (dimensions  in  meters)  for  a  hot  air  system  with 
pebble  bed  storage.   Researchers  generally  quote  suggested 
ratios  of  collector  size  per  storage  volume  (or  per  person 
for  strictly  domestic  water  heating  applications)  when 
speaking  of  water-type  systems.  When  doing  so,  the 
validity  of  the  ratio  is  predicated  on  specifying  the 
location  and  percentage  to  be  met  of  the  total  demand. 
It  is  for  this  reason,  with  a  similar  claim  of  validity, 
that  the  ratio  above  is  presented. 
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APPENDIX 

The  flow  chart  included  here  illustrates  the  simulation 
of  the  solar  heating  system  used  in  this  study. 

The  numbers  appearing  in  the  chart  are  keyed  to  the 
various  "go  to"  statements  used  in  the  computer  program  and 
are  included  in  order  to  guide  the  reader  through  the 
computer  listing  appended. 

The  terms  used  herein  are  defined  in  the  SIMULATION 
NOMENCLATURE  section  of  the  appendix. 

Two  pages  of  output  are  presented  to  illustrate  their 
form. 
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TABLE   A.l 

Sample  of  Yearly  Simulation  Output  for  an 
Average  Day  in  Each  Month 


Collector  Area:   5m  x  12m 

2 
Storage  Volume:   6.8m  x  2.5m 
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COMPUTER  NOMENCLATURE 
FOR  FIGURE  A. 2 


AREAC 

DEGDAY 

AT 

At 
stor 

AQPB 

AQRB 
DLTIME 

DLTIMH 
IFLAG 


QADINT 

QADTOT 

QAUX 
QAXQRM 

QI 
QLOSSH 


surface  area  of  collector 

heat  load  of  dwelling;  (W-HR/°C-HR) 

difference  in  temperature  between  storage  bed 
surroundings  and  average  temperature  in  bed 

change  in  temperature  of  storage  due  to 
heating  bed  with  the  collector 

heat  transferred  from  storage  while  storage 
is  idle  (convective  losses) 

heat  transeferred  to  storage  from  collector 

time,  in  seconds,  remaining  in  any  given  hour 
for  performing  any  of  the  several  modes  of 
heat  transfer  possible  in  this  simulation 

one  hour  in  seconds 

when  equal  to  1,  indicates  the  collector 
is  presently  heating  storage  after  having 
heated  the  dwelling;  used  as  an  internal 
logic  signal  only 

heat  transferred  in  an  interval  (either  60  sec 
or  5  sec)  when  storage  used  to  heat  dwelling 

heat  transferred  while  heating  room  from 
storage 

auxiliary  heat  required  by  dwelling 

monthly  ratio  of  auxiliary  heat  supplied  to 
room  needs 

solar  radiation  incident  upon  collector 


Total  amount  of  heat  lost  from  storage  to 
room  whiie  storage  being  heated  by  collector 
(QLOSS  refers  to  the  amount  lost  in  either 
60  sec  or  5  sec) 
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QRATIO 

QRMFMC 
QRMFMS 
QROOM 
QROOM1 

QSLRQR 

QSTFMC 
QSUBLS 
QU 

QUQI 

QUQRM 

QXCES 


QXCSMT 
SI 

SURFAR 

TAMB 

TAVCHK 

TREF 

TSAVG 


ratio  of  energy  supplied  from  storage  in 
small  interval  to  the  room  needs  in  that 
interval 

heat  transferred  to  the  dwelling  from  collector 

heat  transferred  to  the  dwelling  from  storage 

hourly  room  energy  requirement 

hourly  room  energy  requirement  accounting 
for  the  uncontrolled  storage  loss 

monthly  ratio  of  solar  energy  supplied  to 
room  needs 

heat  transferred  to  storage  from  collector 

uncontrolled  convective  losses  from  storage 

hourly  thermal  energy  extracted  from  the 
collector  by  the  system  =  QRMFMC  +  QSTFMC 

ratio  of  QU  to  QI 

ratio  of  QU  to  QROOM 

the  hourly  amount  of  energy  supplied  to  the 
residence  for  which  there  is  no  need 
(convective  losses  or  "overshoot"  from 
storage) 

monthly  total  of  QXCES 

flux  of  solar  radiation  incident  upon  the 
collector  surface  area 

surface  area  of  storage  bed  through  which 
convective  loss  takes  place 

ambient  temperature 

maximum  average  storage  temperature  allowed 

average  value  of  allowable  room  temperature 
range 

average  instantaneous  temperature  in  storage  bed 
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LISTING  OF 

COMPUTER 
SIMULATION 

TABLE  A. 3 
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